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For a model glass-former we demonstrate via computer simulations how macroscopic dynamic 
quantities can be inferred from a PEL analysis. The essential step is to consider whole superstruc- 
tures of many PEL minima, called metabasins, rather than single minima. We show that two types 
of metabasins exist: some allowing for quasi- free motion on the PEL (liquid-like), the others acting 
as traps (solid-like). The activated, multi-step escapes from the latter metabasins are found to 
dictate the slowing down of dynamics upon cooling over a much broader temperature range than is 
currently assumed. 



I The theoretical understanding of supercooled liquids 
and glasses is one of the important unsolved problems of 

' statistical physics. Since the pioneering work of Gold- 
stein I*!, it has been realized that the potential energy 
landscape (PEL) viewpoint is useful for characterizing 
supercooled liquids and glasses. To this end, one consid- 
ers the high-dimensional vector of all particle coordinates 

] as a point moving on the surface of the total potential 
energy At sufficiently low temperatures, the sys- 

tem resides near the local minima of the high-dimensional 
PEL. It has turned out that the PEL description of ther- 
modynamics starts to work when cooling below approx- 
imately T — 2r„, where T„ is the critical temperature 
of mode-coupling theory Q . In this temperature regime 
the statistical properties of PEL minima fully determine 
all thermodynamic properties 0, H 1^ .Ei .10] . 

What is still lacking in the literature is a full quantita- 
tive understanding of dynamics, i.e. of the slowing down 
of molecular motion upon cooling. In this connection, 

' the single-particle diffusion coefficient, D{T), is a typical 
macroscopic quantity of interest. In general, knowledge 
of thermodynamic properties is not sufficient to predict 
D(T). However, the validity of the Adam-Gibbs rela- 
tion [ly , connecting D{T) to the configurational entropy, 

\ suggests the presence of a strong link between thermo- 
dynamics and dynamics [T^ ll^ ll^ . For a model glass- 
former we have demonstrated via computer simulations 
that D{T) can indeed be inferred from a PEL analysis 
in a quantitative way. The essential step, however, is 
to consider whole superstructures of many PEL minima, 
called metabasins .2i,^il5j, J^, JT, 18], rather than single 
minima. Metabasins are reminiscent of protein-folding 
funnels (19| or related structures in small clusters pol ]. 

' The goal of this paper is to derive physical conse- 
quences of this mapping between PEL and dynamics and 
thus to obtain a coherent picture of the glass transition. 
In particular we will dwell on the interpretation of the 
activation energies for fragile systems in the Arrhcnius 
dia gram , the concept of liquid-like and solid-like behav- 
ior [23, and the question of a possible crossover at Tc, 
as implied by the mode-coupling temperature and the 
analogy to p-spin models. 



In what follows we analyze a binary mixture of 
Lennard- Jones particles (BML J) . All details of the BML J 
have been reported in [T3. This system is designed to 
model a Ni8oP20 mixture. Langevin molecular dynamics 
have been employed with a simulation box of = 65 
particles at density p = 1.2a^f. Throughout the pa- 
per, all quantities are given in reduced units, which, 
in the case of NigoP2o, correspond to a,^f = 2.218 A, 
e„f = 7.762 kJ/mol, t„f = 1.323 fs, r„f = 934 K, and 
Djof = fL/^rof = 0.372 cm^/s. For the temperature 
regime of our analysis, a system size of A^ = 65 is prov- 
ably sufficient to reproduce the dynamic properties of 
a macroscopic system We have checked, e.g., 

that D{T) is the same for A^ = 130 within 10% for 
T > = 0.45T,„f. The minima are obtained by a steep- 
est descent minimization. As usual we perform regular 
quenches of the system to monitor the energies of the cor- 
responding energy minima. In addition, we use temporal 
interval bisectioning to resolve the elementary transitions 
between minima |l7|. 

The definition of metabasins (MBs) is motivated by 
the observation that during the molecular dynamics in 
configuration space the system performs several back- 
and-forth jumps between adjacent minima until finally 
this region is left. The description of the transport may 
be simplified if these minima are regarded as a single su- 
perstructure, i.e. a MB H^. Thus, the MBs correspond 
to an appropriate tiling of configuration space. Although 
irrelevant for the understanding of thermodynamic prop- 
erties, they are of outstanding importance for dynamics. 
As shown in [3] it is possible to find a strict definition 
of MBs which can be also used in practice. Similarly 
to the previous consideration of energy minima the time 
evolution of the system may be regarded as a continuous 
sequence of MB visits with individual residence (waiting) 
times r. 

We briefiy summarize the main results of our analysis 
which will be important throughout the paper. They all 
hold for T < 2Tc- (i) Analyzing the escape characteristics 
from MBs by repeated simulations from the same MBs 
at different temperatures T it turns out that the 
mean residence time in MBs of energy e is given to a 
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good approximation by 

(r(e,r)} «To(e)exp(i?(e)/fc3r). 



(1) 



The Arrhenius form of (T(e, T)) suggests that the escapes 
from MBs are thermally activated, with a barrier height 
of E(e), introduced here as a fitting parameter. The 
prefactor To(e) turns out to be basically independent of 
e. (We define e as the energy of the lowest minimum 
within a MB.) (ii) E{e) is directly related to the PEL 
barriers in the high-dimensional surrounding of the MBs. 
By a detailed discussion of this issue, incorporating the 
funnel-like nature of MBs, this was verified jl^]. Thus, 
by analyzing the local topology of MBs, E{e) can be pre- 
dicted. We may interpret E{e) as the depth of a MB 
of energy e. The deeper the MB in the PEL, i.e. the 
lower e, the higher the activation energy. In this analy- 
sis one has to take into account that the escape from a 
MB is a multi-step process, i.e. generally comprises hops 
between several minima, (iii) A crucial quantity is the 
average residence time {t(T)). It is defined as the average 
over all MBs encountered at a specific temperature, i.e. 
{t{T)) = J dTTip{T,T). By (p{T,T) and (^(e,T), respec- 
tively, we denote the distribution of waiting times and en- 
ergies of visited MBs. (iv) The population of MBs with 
energy e is given by p{e,T) = {T{e,T))^{e,T) / {t{T)) . 
Numerically, it is indistinguishable from the population 
of minima, normally studied in this field and is purely 
gaussian in the accessible energy range Using EqQ] 
the inverse average waiting time can be written as 

^o/(r(T))« j depie,T)exp{-Eie)/k^T) 

= eM-E.n{T)/k^T), (2) 

where the energy dependence of tq (e) has been neglected. 
E^h{T) can be interpreted as the typical barrier height 
encountered at temperature T. (v) A simple relation 
exists between (t(T)) and D{T) via 



D{T) 



6N{t{T))'' 



(3) 



with a temperature-independent effective jump length 
a « 1.0cr„f Thus, the temperature dependence of D(T) 
is exclusively determined by the average waiting time, 
whereas all spatial aspects of diffusion are temperature 
independent. Such a simple relation does not hold on 
the level of single minima. Actually, as a side effect it 
turns out that the dynamics can be basically described 
as a random- walk between MBs, thus suggesting a non- 
topographic view. With the definition Dq = o^/GNtq we 
can thus write 



D{T) « Doexp{~E^,,{T)/k^T). 



(4) 



Note that the ingredients of the constant Dq, i.e. a and 
To, have been obtained from the simulations mentioned 
above. In the units of our simulation we have Dq = 
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FIG. 1: Estimating the diffusion coefficient from the depths 
of MBs and their distribution, (a) The barrier height -B(e) 
as a function of MB energy, and the population of MB en- 
ergies, p{e,T), for emin — — SOBcrof at different temperatures. 
p{e,T) is purely gaussian [2^. (b) Comparison of _Doat(r) 
with D{T). The three curves correspond to different lower 
PEL cutoffs e„i„/e„t = -302, -304, and -306, from bottom 
to top. Different extrapolations of E{e) to e < — 302eref (see 
(a)) lead to nearly identical _Do3t(r)'s (no visible difference 
between the curves corresponding to (i) and (ii) in (a)). We 
set Tg ^ 0.3T„f 0. 



o'^/QNtq « 1.3 • lO^^Acf If the BMLJ is mapped on a 
Ni8oP20 alloy this corresponds to Dq k, 4.8 x 10"^ cm^/s. 

Here we discuss the relevant physical implications of 
Eas lll4l In Fig^a) we show the barrier height E{e) as 
a function of MB energy, and the population of MB en- 
ergies, p{e,T), for different temperatures. From these 
quantities we obtain an estimate for the diffusion coeffi- 
cient, denoted D^^t{T). A comparison of the estimated 
with the true diffusion coefficient D{T) (directly from 
molecular dynamics) is shown in (b). The very good 
agreement at T„ < T < 2T„ reflects the consistency of 
our approach and shows that we can indeed understand 
the macroscopic dynamics from knowledge of the ther- 
modynamics and the local barriers. 

From the simulations, we know p(e, T) and E(e) for 
e > — 302e,„f. Both quantities do not display finite-size 
effects in this energy range Due to the Boltzmann 
weighting in Eq[21one may hope that the available infor- 
mation is sufficient to reliably estimate D^^^iT) also for 
T < Tc- To check this hypothesis, we have calculated 
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D^st{T) thereby using extreme candidates for possible e- 
extrapolations (e < — 302e„f); see Fig^ It turns out that 
-DcBt {T) is very insensitive to -B(e) . For p(e; T) we have ex- 
tended the gaussian description but have set p{e; T) — 
below different cutoff-energies e„i„. (A lower bound for 
e„i„ is given by the condition that the number density 
of minima is of the order of l/e„t @- This leads to 
-306e,rf < £„„„ < -302e„f J^].) As shown in Fig[Tl(b), 
the resulting uncertainty in the estimation of D^^-i(T) is 
small up to T > 0.35T^^f and becomes significant for lower 
temperatures. Since a cutoff at — 302e„f is probably far 
too drastic (see for strong indications), however, the 
slowing down can be thus predicted semi-quantitatively 
down to T w Tg, a region which very likely will never be 
accessible by molecular dynamics simulations. 

Moreover, Eq^ implies that the activation energy 
which is obtained from connecting (l/T = 0/T„f, Z) = 
-Do) and {l/T, D{T)) in an Arrhenius plot by a straight 
line has a simple interpretation: it is the typical bar- 
rier height the system experiences at a given tempera- 
ture j27l |. This non-trivial interpretation has very re- 
cently been confirmed by hyperquench experiments |28j| . 
where the typical barrier height was measured by probing 
the specific heat during the reheating process. 

To obtain a more detailed picture of MBs we have de- 
termined the distribution of waiting times, ip{T,T), as 
shown in Fig|21[a). By visual inspection, one can distin- 
guish two different time regimes (separated by the ver- 
tical line at r* k, SOOOt^^f). In terms of energetics one 
expects that the long r's arise from deep traps where the 
system is caught for a long time. In the opposite limit 
one may imagine that there are quite shallow MBs which 
do not strongly confine the system so that it will mainly 
stay close to the high-dimensional boundary of these MBs 
(see the sketches in (a)). Thus, short waiting times cor- 
respond to just scratching the MBs. This expectation 
can be verified by analyzing the trajectory of the system 
during the MB visits; an example of such a computa- 
tion is given in Fig^b). There we plot (cosa(r)) over 
MBs of lifetime r, where a is the angle between the entry 
point, the lowest minimum of the MB and the exit point 
of the trajectory (see the sketches in (b)). Again, two 
time regimes with a temperature-independent crossover 
time T* can be identified. Short visits to MBs lead to 
small values of a, meaning that the system indeed merely 
scratches these MBs. For r > r* the value of (cos a(r)) 
reaches a limiting value of ca. 0.2 , indicating that entry 
and exit points are largely uncorrelated. This should be 
the case after a long equilibration inside a MB with many 
possible exits. Due to the difference in stability, we call 
the MBs with r > r* solid-like, the other MBs liquid-like. 
This notation has been borrowed from two-state models 
where these two types of configurations have been pos- 
tulated [2ll |. 

Clearly, the molecular slowing down upon cooling is 
caused by the enhancement of solid-like configurations. 
This can be quantified in two different ways. Firstly, 
one may determine the fraction of solid-like configura- 




waiting (residence) time in MBs 



FIG. 2; Different exploration of MBs, as evidenced by the 
waiting-time distributions and the correlations between MB- 
entry and exit points, (a) The distribution ipij, T) of wait- 
ing times at different temperatures, (b) Average value of 
(cos Q(r)) over MBs of lifetime r (see text and the sketches 
in (b), x{t) symbolizes the high-dimensional trajectory of the 
system) . 



tions the system encounters, i.e. ip^ai{T) = dTip{T, T), 
where ip^^i < 0.5 implies that more liquid-like than solid- 
like configurations are visited. Secondly, we can spec- 
ify the fraction of time spent in solid-like configurations, 
which can be expressed as p^atiT) = dTp{T,T) with 
p{t,T) = t(p{t,T)/{t{T)). > 0.5^implies that the 
system is mostly residing in solid-like configurations. In 
Figl^l we show the temperature dependence of ip^^i and 
Paoi- Three different temperature regimes can be distin- 
guished (see the sketch in FigOJ: For T > 2T, both 
quantities are smaller than 0.5. Thus, the system be- 
haves liquid-like. Interestingly, this temperature regime 
(defined by dynamics) is exactly the temperature regime 
for which the minima no longer influence the thermo- 
dynamic properties of the system Q. Below 2T^, p^^i is 
larger than 0.5. Thus, the system mainly resides in solid- 
like MBs. Finally, below a temperature near T^, also (p^^i 
exceeds 0.5, i.e. we have a trap-to-trap motion. This 
crossover, however, is very gradual. 

With Eq[21we can also analyze the question in which 
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FIG. 3: The temperature dependence of psoi (the fraction of 
time spent in soUd-like configurations), (the fraction of 
solid-like configurations encountered during the time evolu- 
tion), and Pact (a measure for the contribution of activated 
processes to the diffusion coefficient; see text). In the lower 
part we depict schematic plots of the scenarios in the three 
temperature regimes. The squares symbolize the different 
MBs. 

temperature regime the dynamics are dominated by ac- 
tivated processes. We call a process activated if the ac- 
tivation energy E{e) is larger than 5ksT. To this end, 
we consider the fraction of (t(T)) which is made up by 



activated hops, i.e. 

(r(r))-i / de(T(e, T))^(e, T) ^ p^T), 

JE(e)>5kuT 

as shown in Fig|3| The crossover temperature for which 
Pact = 0.5 is close to 1.5T„. Thus, already signifi- 
cantly above T^, the dynamics are dominated by acti- 
vated processes. Please note that the MBs contribut- 
ing to the above i nteg ral are solid-like since they fulfill 
(T(e,T)) > T* (see This is why Pact(r) < p,„,(T). 

The present results seem to be at variance with the 
current understanding that activated transitions between 
PEL minima set in at T^, whereas above T„ they are irrel- 
evant. The latter scenario is motivated by the properties 
of p-spin models and has been backed by the ob- 
servations that the number of free directions (obtained 
via instantaneous-normal- mode analyses [scl l3l|). is di- 
rectly related to the diffusion coefficient and that above 
the system is close to transition states rather than 
close to minima [s^, H^i 13- In our earlier publica- 
tion , we have already pointed out the technical short- 
comings of the latter kind of investigations. More impor- 
tantly, though, the above-cited works do not incorporate 
the MB structure of configuration space. In doing so, one 
has focused on the vast majority of mira-MB transitions, 
which, however, are irrelevant for relaxation. 

Finally we note that the results of our work with re- 
spect to the crossover to activated behavior as well as the 
non-topographic nature of inter-MB dynamics is compat- 
ible with the implications obtained from a recent analysis 
of spin- facilitated models of the glass transition js^j ]. 

We thank C.A. Angell, L. Berthier, J. P. Garrahan, and 
H.W. Spiess for helpful discussions. 
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